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Cygnus X-l was the first X-ray source widely accepted to be a black hole candidate 
and remains among the most studied astronomical objects in its class. The detection of 

Q ■ non-thermal radio, hard X-rays and gamma rays reveals the fact that this kind of objects 

5-H ' are capable of accelerating particles up to very high energies. 

r/^ ' In order to explain the electromagnetic emission from Cygnus X-l in the low-hard 

C^ , state we present a model of a black hole corona with both relativistic lepton and hadron 

content. We characterize the corona as a two- temperature hot plasma plus a mixed non- 
thermal population in which energetic particles interact with magnetic, photon and mat- 
ter fields. Our calculations include the radiation emitted by secondary particles (pions, 
muons and electron/positron pairs). Finally, we take into account the effects of photon 
yK-* ■ absorption. We compare the results obtained from our model with data of Cygnus X-l 

U^ ' obtained by the COMPTEL instrument. 

m 

l/) ' Keywords: black hole; corona; gamma-rays. 

in 

f~3 1. Basic scenario 

The low-hard state of accreting black holes is characterized by the presence of a hot 
corona around the compact object. Figure [1] shows a scheme of the main components 
of the system. For this geometry, we assume a spherical corona with a radius R c and 
an accretion disk that penetrates the corona up to a radius R^ < R c . We suppose 
that the corona is homogeneous and in steady state. 

We assume that the luminosity of the corona is 1 % of the Eddington luminosity 
of the black hole. For a corona characterized by the parameters of Table[TJ this yields 
L c = 1.3 x 10 37 erg s" 1 . The corona is composed of a two-temperature plasma, with 
an electron temperature T e = 10 9 K and an ion temperature Tj = 10 12 K (e.g. Ref. 
jTJ) . This is a consequence of electron cooling being more efficient than proton. 

In the low-hard state we can assume the system is in a steady state. We then 
consider equipartition of energy between the different components of the system 
and obtain the values of the relevant physical parameters. 

The X-rays emission of the corona is characterized by a power law with an 
exponential cut-off at high energies, 

(A ph E- a E<E C 

n M^) - | AphE -<* e -E/E c E >E C . 



June 1, 2010 0:16 WSPC/INSTRUCTION FILE Vieyro 



2 F.L. Vieyro, G.E. Romero & G.S. Vila 




Accretion 
disk 



Fig. 1. Schematic representation of the components of the system. (Not to scale) 

In accordance with what is observed in Cygnus X-l (e.g. Ref. [5]), we adopt a = 1.6 
and E c = 150 keV. 

Table Q] shows the values of the parameters in our model. 

Table 1. Model parameters. 



Parameter 


Value 


Mbh : black hole mass [Mq] 


lOf 1 ) 


R c '- corona radius [cm] 


5.17 x 10 7(1 ' 2) 


T e : electron temperature [K] 


10 9 


T;: ion temperature [K] 


10 12 


E c : X-ray spectrum cut-off [keV] 


150 


a: X-ray spectrum power-law index 


1.6 


r\: acceleration efficiency 


10- 2 


B c : magnetic field [G] 


2 x 10 7 


rii,ne: plasma density [cm -3 ] 


6.2 x 10 13 


a: hadron-to-lepton energy ratio 


1 - 100 


kT: disk characteristic temperature [keV] 


0.1 


v. advection velocity [c] 


0.1 



C 1 ) Typical value for Cygnus X-l (Ref.|2). 
< 2 ) 35i? G . 



2. Particle injection 

The injection function for relativistic protons and electrons is a power-law in the 
energy of the particles Q(E) — QqE~ y e~ ~ E '' ' £m " , which is the consequence of dif- 
fusive particle acceleration by shock waves in the coronal plasma. Typical spectral 
indices are T ~ 2.2 (Ref. [3}. The normalization constant Qq can be obtained from 
the total power injected in relativistic protons and electrons, L re i = L p + L e . This 
power L rc i was assumed to be a fraction of the luminosity of the corona, L rc \ = kL c , 
with k = 10 -2 . The way in which energy is divided between hadrons and leptons 
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is unknown, but different scenarios can be taken into account by setting L p = aL e . 
We consider models with an injection dominated by protons in which a = 100 and 
models with a = 1. 

3. Particle losses and acceleration 

The injected particles will lose their energy by different radiative processes. There 
are three relevant processes of interaction of relativistic electrons with magnetic, 
matter and photon fields in this scenario: synchrotron radiation, inverse Compton 
scattering, and relativistic Bremsstrahlung. For protons there are also three rele- 
vant processes: synchrotron radiation, proton-proton inelastic collisions and photo- 
hadronic interactions. These interactions produce pions, which then decay producing 
muons. We also consider the energy loss of these particles; the processes for charged 
pions are the same as for protons and for muons the same as for electrons. 

We consider two mechanisms of particle escape from the corona: advection and 
diffusion. In the case of advection, particles fall onto the compact object at a mean 
radial velocity v = 0.1c (Ref. 3). Therefore the advection rate is 

a = £• (i) 

In the case of diffusion, we consider that the corona is static and diffusion of relativis- 
tic particles occurs in the Bohm regime. The diffusion coefficient is D(E) = r g c/3, 
where r g = E/(eB) is the giroradius of the particle. The diffusion rate is 

*diff- R 2 ■ W 

Figure [2] shows the cooling rates for different processes of energy loss, along 
with the acceleration and escape rates, for each type of particle considered. Under 
the physical conditions previously described, the main channel of energy loss for 
electrons is synchrotron radiation. For protons, both pp and pry interactions are 
relevant. However, in the model with advection, most protons fall into the black 
hole before radiating their energy. For pions, the main channel of energy loss is 7r7 
interaction, but an important fraction of pions decay before cooling. On the contrary, 
muons with energies above 10 11 eV cool mostly by inverse Compton scattering. 

The maximum energy that a relativistic particle can attain depends on the 
acceleration mechanism and the different processes of energy loss. The acceleration 
rate t^r c for a particle of energy £ in a magnetic field B is given by 

C = ^, (3) 

where r\ < 1 is a parameter that characterizes the efficiency of the acceleration. We 
have fixed r\ = 10~ 2 , which describes an efficient acceleration. 
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Fig. 2. Radiative losses in a corona characterized by the parameters of Table [T] 



4. Spectral energy distributions (SEDs) 

In order to obtain the spectral energy distributions produced by the different ra- 
diative processes, we solve the transport equation for each kind of particle. We cal- 
culate the SEDs of different processes following Ref. ((5|) as standard reference. We 
also calculate the radiation emitted by secondary pairs, which are injected mainly 
by photon-photon annihilations. 

Figure [3] shows all contributions to the total luminosity for different advec- 
tion regimes and for different values of the parameter a. Luminosities produced by 
hadrons and muons are higher in models with diffusion. This is because in models 
with advection an important fraction of protons and pions are swallowed by the 
black hole, while with diffusion these particles are able to lose their energy before 
falling into the compact object or escaping. 

Finally, we take into account the effects of photon-photon attenuation. The ab- 
sorption can be quantified through the opacity r. If the initial gamma-ray luminosity 
is L^(Ej), the attenuated luminosity L 1 (E~ f ) after the photon travels a distance I 



-<y\,— '7 

is L 1 (E, 



LJE~ 



,-r(l,E^ 
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Fig. 3. Spectral energy distributions obtained for different set of parameters. 



Figure S] shows the attenuated luminosity for each set of parameters, the spec- 
trum of Cygnus X-l as observed by COMPTEL (Ref. 6), and the sensitivity of 
different instruments. Notice that in addition of accounting for the MeV tail ob- 
served by COMPTEL, a high energy bump with luminosities ~ 10 33 erg s _1 is 
expected at E 1 > 1 TeV. 

The existence of the high-energy tail in the spectrum detected by COMPTEL 
and reproduced by our model was recently confirmed by observations with INTE- 
GRAL satellite (Ref. [7]). 



5. Conclusions 

On one hand, the SED predicted by the model with an injection dominated by 
protons and diffusion reproduces the observations of Cygnus X-l. On the other 
hand, our model predicts a high energy emission, above 1 TeV, that might be 
detectable by instruments like MAGIC II or CTA in the near future. 
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Fig. 4. Spectral energy distributions. Observational data from COMPTEL of Cygnus X-1. 
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